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Abstract 

We present estimation of proton-air inelastic interaction cross-section obtained for different energy values 
are. Results are compared with different hadron interaction models. 
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I. INTRODUCTION 



For more than 40 years the Yakutsk array continuously records ukra-high energy cosmic rays 
(UHECR) by measuring electromagnetic, muonic and Cherenkov components of extensive air 



showers (EAS) UJ]. Throug 



1 these measurements a lateral distribution of charged particles and 



muons with cthr. > 1 GeV [[21, sj] was studied and longitudinal development of air showers was 
reconstructed by measuring the lateral distribution of Cherenkov light emitted in EAS 

Estimation of the depth of maximum shower development (Xmax) in individual showers with 
energy 10^^ — 10^^ eV at the Yakutsk array was performed using readings from integral and differ- 
ential (tracking) Cherenkov light detectors. Integral detectors provide observation of Cherenkov 
photons lateral distribution in wide distance range, while differential detectors record time-base 
of the Cherenkov light pulse at different core distances. It follows from model calculations that 
characteristics such as P = logio <5(200)/(5(550) (a relation between the density of Cherenkov 
light fluxes at 200 and 550 m from shower axis) and ri/2 (the width of the Cherenkov pulse at the 
half-height) are sensitive to the depth of maximum shower development. Besides, the P parame- 
ter depends weakly on conditions of Cherenkov light distribution in the atmosphere (~ 8 %), but 
depends strongly on the accuracy of shower geometry reconstruction (~ 13 %). The ri/2 is mea- 
sured by equipment that provides high precision of the Cherenkov pulse shape reproduction and 
thus assures better accuracy in a;niax measurement (~ 25 g/cm^^). The methodology of Cherenkov 
light measurement and analysis is presented in a greater detail in papers [3, 6]. 



n. METHODOLOGY OF THE ANALYSIS 

To extract information about the parameters of interaction between the primary particle and 
air nuclei it is necessary to measure those characteristics of EAS which are morphologically con- 
nected to the start of nuclear cascade in air. One can test the sensitivity of EAS characteristics 
to cross-section of inelastic interaction with Monte-Carlo method to simulate artificial showers 
using some model of hadronic interactions. Such calculations were performed in multiple works. 
Further we will consider the methodology and algorithm for estimation of inelastic interaction 
cross-section involving the data on the depth of maximum EAS development (a;niax)- 

Total cross-section of secondary particles production in EAS is one of fundamental values that 
characterize hadronic interactions. For showers initiated by primary protons of some energy it is 

2 



TABLE I. Cross-sections of p-air inelastic interaction 



£^0, eV Aexp. ki k2 Aphys., g/cm"^ o-PJ^°^;, mb 

3.3 X 10^^ 110.3 ± 4.1 1.4 ± 0.3 1.25 63.0 ± 5.2 385 ± 23 

7.2 X 10^5 107.4 ± 4.1 1.4 ± 0.3 1.25 61.4 ± 5.2 396 ± 23 

2.1 X 10^6 102.6 ± 5.2 1.4 ± 0.3 1.25 58.7 ± 5.2 413 ± 34 

possible to obtain the distribution of the depth of initial interaction (xi): 



dp 1 ^ 

e , (1) 



dXi Ap.air 

where Ap.air is a mean free path of a proton before interaction in air. Mean depth of first interaction 
with respect to its fluctuations is directly connected to the cross-section of inelastic interaction 
between proton and air nucleus: 

O^p-air = -T (^) 

^p-air 

Here, (rriair) is mean air "mass", which is ~ 14.5mp ~ 24253.01 mb/gcm^ [7]. The relations 
© and © were used to estimate the inelastic interaction cross-section by experimental data ob- 
tained at the Yakutsk array. During the measurement of the Ap.air value, a mass composition of 
primary cosmic rays (PGR) was taken into account (contribution to generation of Xmax distribution 
at fixed energy of other nuclei) by utilizing the points shifted by one a(xmax) from the maximum 
of histogram (see figlD- 



III. EXPERIMENTAL DATA 

In present work we employ measured Xmax distribution at fixed energy values 3.3 x 10^^ eV, 
7.2 X 10^^ eV and 2.1 x 10^^ eV which were obtained from Cherenkov light data recorded by 
the small Cherenkov setup during more than ten years of continuous observation (see figUl)- Final 
results of the analysis are presented in table H] and fig HI 

FiglD shows a compilation of the (Jp.air(-Eo) data obtained in different experiments including 
the Yakutsk array in the energy range 10^'' — 10^^ eV (see Tien Shan data for 1 PeV [9]). In 
this work, using similar algorithm and performing the simulation of x^ax measuring process with 
Monte Carlo method (with respect to hardware uncertainties), we obtained the ctP™'' for energies 
10^^ — 10^^ eV. New results are shown on figj3l On the same figure results of calculations with 



various hadron interaction models are shown 



]. Extrapolation of accelerator data leads to a large 
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FIG. 1. Fluctuations of Xmax at fixed energy and tlie teclinique of cjp.air measurement by Clierenkov light 
data adopted at tlie Yakutsk array. 983 events, Eq = 2.1 x 10^'' eV, Agxp = 102.6 ± 5.2 g/cm"^, x^ax = 
572 ± 6.6 g/cm-2, ^(xmax) = 56 ± 3 g/cm^^. 



uncertainty at energy above 10^^ eV. To a greater extent this is tied to some discrepancy between 
models in description of more precise accelerator data and different approaches in description 
of generation of secondary particles in EAS. As it is seen on fig. [3l cosmic ray data together 
with accelerator data allow to introduce some limitations on applicability of particular hadron 
interaction model for description of ultra-high energy air shower development. 



IV. CONCLUSION 

The Yakutsk EAS array is a medium sized array compared to compact and large arrays. It ef- 
fectively controls the energy range from 10^^ eV to 5 x 10^^ eV. In several decades of continuous 
observations a large amount of experimental data was accumulated. Both irregularities in cosmic 
ray energy spectrum are manifested in these data and the mass composition of primary particles 
initiating ultra-high energy air showers has been estimated using one technique. Such simulta- 
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FIG. 2. Initial Xmax distributions, from which the Agxp. was derived. 

neous analysis of the spectrum and mass composition allows one to test various hypotheses on 
cosmic rays origin and their propagation in the Universe. To increase the precision of this analysis 
it would take a choice of a single hadron interaction model from the group of models presented to- 
day. Studying the energy dependence of inelastic interaction cross-section just helps make such a 
choice. For example, the choice of the model for EAS development would allow experimenters to 
improve the precision in energy estimation and obtain the model-independent estimation of cosmic 
ray mass composition in ultra-high energy domain, more reliable than those adopted nowadays. 

The work is supported in part by SB RAS (integral project "Modernization of the Yakutsk 
array"), RFBR (grant#l 1-02-00158) and the Russian Ministry of Education and Science (contract 
#02.740.11.0248). 
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FIG. 3. A comparison of current results with the other experiments and with some hadron interaction 
models. The original figure is taken from the work 171]. 
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